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Abstract 

Previously Chicone, Latushkin and Montgomery-Smith [Comm. Math. Phys. 173,(1995)] 

have investigated the spectrum of the dynamo operator for an ideally conducting fluid. More 
recently. Tang and Boozer [Phys. Plasmas (2000)], have investigated the anisotropies in 
magnetic field dynamo evolution, from finite-time, Lyapunov exponents, giving rise to a Rie- 
mann metric tensor, in the Alfven twist in magnetic flux tubes (MFTs). In this paper one 
investigate the role of Perelman Ricci flows constraints in twisted magnetic flux tubes, where 
the Lyapunov eigenvalue spectra for the Ricci tensor associated with the Ricci flow equation 
in MFTs leads to a finite-time Lyapunov exponential stretching along the toroidal direction 
of the tube and a contraction along the radial direction of the tube. It is shown that in the 
case of MFTs, the sectional Ricci curvature of the flow, is negative as happens in geodesic 
flows of Anosov type. Ricci flows constraints in MFTs substitute the Thiffeault and Boozer 
[Chaos(2001)] have vanishing of Riemann curvature constraint on the Lyapunov exponential 
stretching of chaotic flows. Gauss curvature of the twisted MFT is also computed and the 
contraints on a negative Gauss curvature are obtained. 

PACS numbers:02.40.Hw:differential geometries. 91. 25. Cw-dynamo theories. 
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I Introduction 



After the Perelman's seminal paper [1], on an explanation of the Poincare conjecture, 
many applications of the Ricci flow equations to several areas of physics have appeared 
in the literature [2, 3]. More recently Cao [3] has computed the first eigenvalues, of "heat 
flow" operator —A + where R is the Ricci curvature scalar and A := is the Lapla- 
cian 3D operator, under Ricci flow. He found that the eigenvalues were non-decreasing as 
happens in general with Lyapunov exponents [4]. In this paper one investigates a similar 
problem in the context of the time evolution operator under Ricci flow in the backyard 
of the twisted MFT. To be able to determine the Lyapunov eigenvalue spectra in twisted 
MFTs, which is a fundamental problem which helps one to determine the stability of the 
flows inside flux tubes such as tokamaks or stellarators in plasma physics [5] or even in 
the context of solar or other stellar plasma loop [6], one computes the eigenvalue spec- 
tra of the Ricci tensor under the Ricci flow. This is an important problem for plasma 
theorists and experimentalists, shall be examined here in the framework of Riemannian 
geometry [7]. Recently, Thiffeault, Tang and Boozer, investigated Riemannian constraints 
on Lyapunov exponents [8], based on the relation between the Riemann metric and the 
finite-time Lyapunov exponential stretching, so fundamental for dynamo action. Note 
that, here, another sort of constraint is investigated. Instead of the vanishing of the Rie- 
mann curvature tensor, called by mathematicians, Riemann-flat space or condition, one 
use the Lyapunov spectra under Ricci flow. Chaotic flows inside the twisted MFTs are 
investigated. Anosov diffeomorpism [9], is an important mathematical tool from the the- 
ory of dynamical systems, that has often been used, in connection with the investigation 
of dynamo flows and maps [10] such as the Arnold's Cat Map [10] on the torus, useful 
in mixing [11] problems in the physics of fluids [11]. One of the main properties of the 
Anosov maps is that they yield Lyapunov exponential of the chaotic exponential stretch- 
ing, which are constant everywhere [12]. This paper is organized as follow: In section 2 
the dynamo maps under Ricci flows in MFTs are investigated with the aid of Lyapunov 
spectra. In section 3, the thin tube perturbations are computed in the negative Ricci 
sectional curvature assumption. Section 4 addresses discussions and conclusions. 
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II Ricci fast magnetic dynamo flows in MFTs 



Let us start this section, by defining the Ricci fiow as: 
Definition 2.1: 

Let us consider a smooth manifold which Ricci tensor Ric obeys the foUowing equation: 

= 2Ric (ILl) 
dt ^ ' 

Here g is the Riemann metric over the manifold M. where in g(t), te[a, h]. By chosing a 
local chart lA on this manifold, the Ricci flow equation may be written as 



= -2i?i,- (11.2) 



dt 

with this equation in hand let us now compute the eigenvalue spectra of the Ricci tensor 
as 

RijX' = AXi (11.3) 

where here (i, j — 1, 2, 3). Substitution of the Ricci fiow equation (IL2) into the eigenvalue 
equation (11.3) one obtains an eigenvalue equation for the metric itself, as 

= -2A^,,x^- (IL4) 

and since the eigenvector is in principle arbitrary, one can reduce this equation to 



= -2A5,, (IL5) 



dt 

which reduces to the solution 

Qij = exp[-2\it\5ij (IL6) 

one notes that the, 6ij is the Kroenecker delta diagonal unity matrix. No Einstein sum 
convention is being used here. By considering the Tang-Boozer relation between the 
metric gij components and the Lyapunov exponents 

Qij = Aieiei + A2e2e2 + Aaeaes (n.7) 
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the following lemma, can be proved: 
lemma 1: 

If Xi is the eigenvalue spectra of the Ricci tensor Ric under Ricci flow equation, the 
Lyapunov spectra is given by the following relations: 

Xi = -7^ < (II.8) 
where Aj are the finite-time Lyapunov numbers. The infinite or true Lyapunov number is 

Here one has used the finite-time Lyapunov exponent given by 

A. = (^) (11.10) 

Note that one of the interesting features of the Ricci flow method is that one may find 
the eigenvalue Lyapunov spectra without computing the Ricci tensor, of the flux tube, 
for example. Actually the twisted flux tube Riemannian line element [13] 

df = dr^ + r^dOn^ + K^{r, s)ds' (ILll) 

can now be used, to compute the Lyapunov exponential stretching of the flow as in 
Priedlander and Vishik [14] "dynamo flow" with the Ricci flow technique above. Here 
twist transformation angle is given by 

e{s) --^Or- j T{s)ds {11.12) 

One another advantage of the method used here is that this allows us to compute the 
finite-time Lyapunov exponential, without the need of recurring to the non-Anosov maps 
[15]. The Lyapunov exponents here are naturally non-Anosov since the exponents are 
non-homogeneous. Here K{r, s) :— (1 — r«;(s, t)cosO), is the stretching in the metric. Let 

us now compute the eigenvalue spectra for the MFTs. Note that the eigenvalue problem, 
can be solved by the 3D matrix 

/2Xgn \ 



MsD 



dtg22 + Xg22 
V dtQss + XgsaJ 



(H.13) 



The eigenvalue equation 

Det[M^j^] = (11.14) 

the following eigenvalue Lyapunov spectra for thick tubes, where K ^ —Korcos9{s) are 

Ai = 0, A2 = 2^, A3 = + coitgeis) (11.15) 

Therefore, since for the existence of dynamo action, at least two of the Lyapunov ex- 
ponents have to have opposite signs [15] in order to obey the stretching (A3 > 0) and 
contracting (A2 < 0), in order that the radial flow be negative, A3 > 0. Therefore from 
the above expressions the following constraint is obtained: 

\ijJxt9Q{s)\ > IA2I = I — I (11.16) 
r 

Thus the ^ < yields a compression on the flux tube which induces the in the tube a 
stretch along the toroidal direction-s, by the stretch-twist and fold dynamo generation 
method of Vainshtein and Zeldovich [16]. In the next section one shall compute the 
relation between twist or vorticity, and the when the sectional curvature of the thin tube 
is negative. Actually this leads us to the following time dependence of the magnetic field 
components, as 

Be « e^^"* = e(v)* (11.17) 
^ e^=* = e^T+'^i'''^*""^)* (11.18) 
where u>i is a constant vorticity inside the dynamo flux tube. 

Ill Magnetic dynamo flows of negative sectional cur- 
vature 

In this section the Ricci sectional curvature [17] in the case of a radial perturbation of a 
thin twisted MFT. This is justified since, as one has seen in the last section, the radial flow 
is fundamental for the existence of non- vanishing Lyapunov exponential stretching, which 
in turn are fundamental for the existence of dynamo action. On the other hand, following 
work by D. Anosov [9], Chicone and Latushkin [18] have previously shown that geodesic 
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flows, which possesses negative Riemannian curvature has a fast dynamo action. Let X 
and Y be vectors laying in tangent manifolds TM. to a Riemannian manifold M. C Af 
where Af is an Euclidean three dimensional space. The Ricci sectional curvature is given 

where R{X, Y)Z is the Riemann curvature given by 

R{X, Y)Z = VxVyZ - VyVxZ - V^x,Y]Z (III.20) 

where 

S{X,Y) := ||X|H|F||2- < X,r >2 (III.21) 
As usual V xY is the Riemannian covariant derivative given by 

VxY = {X.V)Y (III.22) 

Expression [X, Y] is the commutator, where which on the vector frame e;, (/ = 1, 2, 3) in 
R^, as 

X = XkBk (III.23) 

or its dual basis 

X = X^dk (III.24) 
where Einstein summation convention is used. Thus the commutator is written as 

[X,Y] = [X,Yfdk (III.25) 

Thus the Riemann curvature tensor becomes 

R{X, Y)Z = [R^ jkpZ^ X^YP]di (III.26) 

By considering the thin tube approximation K 1 where the gradient is given by 

V = [dr,r-^de^,ds] (III.27) 

Let the radial perturbation in the manifold of flux tube be given by 

X = v\er (in.28) 



where it was considered that the background radial flow vanishes. The vector Y is given 

by 

Y = veee + Vst (III.29) 

where the Frenet vector t is tangent to the magnetic tube axis, and the perturbation 
is chosen along the radial direction, because in general flux tubes initially possesses a 
strongly confined magnetic field along the tube. In the way of computing the Ric, the 
covariant derivative is 

VxY = v\dr[veee + ^;,t] (III.30) 
which, under the approximations 

v\drV0 ^ (III.31) 
v\drVs ^ (III.32) 

The MFT relation 

Vy^xY ^ (III.33) 

Since the term 

V[x,Y]Y = 0{v^) (III.34) 

and one is assuming that the velocities involved in the plasma dynamo fiow are nor small 
neither turbulent, this term can be dropped and the Ric can be expressed as 

R{X, Y)Y ^ VxVyy (III.35) 

This yields 

R{X, Y)Y fa [vs — T{s)~^][v0KTsin9ee — rvosinOt + VsKu] (III. 36) 

Therefore after some algebraic manipulation one obtains the sectional curvature as 

K{X,Y):^^^P^ 

V^rVs 

Where one has considered the approximation that due to dynamo action Vg » vg which 
is usual, for example, in solar physics. In the above computations one has considered the 
fiow as incompressible, V.v = 0, or obeying the equation 

dsVg :— nrrsinOve (III. 38) 



The Gauss curvature is given by 

Ka := (III.39) 

Here g :— detgij and i?i2i2 is the Riemann curvature of the twisted MFT surface, given 
by the hne element 

dl^ = ro^dOR^ + K\s)ds^ (III.40) 

This metric form has been obtained from expression (11.11) by simply consider the Rie- 
mann curved twisted flux tube surface of constant cross-section of radius, r — ro — 
constant. The Riemann curvature i?i2i2 of this Riemannian line element is 

i?i2i2 = -k(s)K(s)cos0 (III.41) 

Here one have considered that the stretching metric coefficient is K[s) = {1 — roK{s)cos6) . 
Thus since g — ro^, substitution of this values into the Kq Gauss curvature above yields 

Ka = (111.42) 

This is a general form of the Gaussian curvature, which can now be applied to the case of 
thin constant Prenet curvature MFTs, by making K 1 and k. — Kq — constant, which 
yields 

Ka = (lllAS) 

Note from this expression that, for it be negative, the cos9 > and the Prenet curvature 
of the magnetic flux tube axis, positive, or they have to have the same sign at all. Just for 

comparison one mention here the fast kinematic dynamo eigenvalue spectrum, in diffusive 
media obtained by Chicone and Latushikin [18], which is 

A, = i [-e(l + k"") + ^e2(l - k2)2 _ 4k)] (HI.44) 

where e is the resistive plasma coefficient of the diffusive fast kinematic dynamo represents 
the growth rate of the magnetized plasma dynamo [19], and In the case of diffusive-free 
(e = 0) ideal plasma one obtains 

Ao = iWA (111.45) 

which is a pure magnetic unity. This is similar to the A3 one has obtained from the 
eigenvalue matrix above that the A3 = IktcjosQ it is also proportional to the Prenet 
curvature. 
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IV Conclusions 



Let us state here for discussion the Chicone-Latushkin fast dynamo in geodesic flow 
theorem [18] as: 

Theorem: If v is the vector field that generates the geodesic flow for a closed two 
dimensional Riemannian manifold M. . then v is a steady solution of Eulers equation on 
M.. In addition if M has constant negative curvature k, then for each magnetic Reynolds 
number Rcm > \/— ^, the corresponding dynamo operator has a positive eigenvalue given 
by Ag above. 

Therefore one may say that in this paper we gu=ive an example of the validity of Chicone- 
Latushkin theorem on geodesic dynamo flows in the case of flux tubes where a containing 
R icci plasma flow. Also in this paper the importance of investigation the Ricci flows as a 
constraint to dynamo flows inside twisted MFTs is stressed, given examples of the dynamo 
action existence in negative sectional curvature. This can be basically done with the help 
of the Lyapunov spectraof the Ricci fast dynamo flows, which are fundamental for the 
exponenial stretching which are in turn so important for dynamo action. An instability 
of the Lyapunov exponential stretching influence on the instability of the Euler equations 
has been discussed by Friedlander and Vishik [14]. Physical implications to the ideas 
discussed here to the dynamo magnetic flux tube discussed and developed by Schuessler 
[?] may appear elsewhere. 



V Acknowledgements 



I am very much in debt to Jean Luc Thiffeault for calling to my attention, many aspect 
of non-Anosov maps and stretching in chaotic flows. I also am deeply indebt to G. 
Paternain , R. Ricca and Dmitry Sokoloff for helpful discussions on the subject of this 
paper. Financial supports from Univcrsidadc do Estado do Rio de Janeiro (UERJ) and 
CNPq (Brazilian Ministry of Science and Technology) are highly appreciated. 



9 



References 

[1] G. Perelman, The entropy formula for the Ricci flow and its geometric applications, 
Los Alamoa Arxives: math/0211599. 

[2] V. Dzhunushaliev, Quantum wormholes as a Ricci flow, Arxiv:0809.0957. 

[3] X. Cao, First eigenvalues of geometric operators in Ricci flows, Arxiv:0710.3947. 

[4] J-L. Thiffeault and A. H. Boozer,Chaos (2001). 

[5] X. Z. Tang and A H Boozer, Phys. Plasmas 11 . 

[6] S.Vainshtein, A. Bykov and I.N. Toptygin, Turbulence, Current Sheets and Shock 
waves in Cosmic Plasmas, Gordon. 

[7] E. Cartan, Riemannian geometry in orthonormal frames, Princeton University Press, 
(2001). A H Boozer, ApJ 394 357 (1992). 

[8] X. Z. Tang and A H Boozer, Chaos 11 (2001). 

[9] D.V. Anosov, Geodesic Flows on Compact Riemannian Manifolds of Nega- 
tive Curvature, (Steklov Mathematical Institute, USSR) vol.90, eds. I. Petrovskii 
and S. Nokol'skii, AMS, Providence, Rhode Island, (1969). 

[10] V. Arnold, Ya B. Zeldovich, A. Ruzmaikin and D.D. Sokoloff, JETP 81 ,n.6, 2052 
(1981). V. Arnold, Appl Math and Mech 36, 236 (1972). V. Arnold and B. Khesin, 
Topological methods in Hydrodynamics (1998) Springer. 

[11] R. Sturman, J. M. Ottino and S. Wiggins, The Mathematical Foundations of 
Mixing. The Linked Twist Map as a Pciradigm in Applications: Micro to 
Macro, Fluids to solids, (2006), Cambridge. Ya B Zeldovich, A. Ruzmaikin and D 
Sokoloff, The Almight Chance (1990) World Scientific. 

[12] S. Childress and A. D. Gilbert, Stretch, Twist, Fold: The fast dynamo , 

Springer, Berlin, New York, (1995). 



10 



[13] R. Ricca, Solar Physics 172 (1997),241. L.C. Garcia dc Andradc, Phys. Plasmas 14, 
(2007). L C Garcia de Andrade, Phys Plasmas 13 (2006). L. C. Garcia de Andrade, 
Non-holonomic dynamo filaments as Arnolds map in Riemannian space, Astronomi- 
cal notes (2008) in press. L. C. Garcia de Andrade, The role of curvature and stretch- 
ing on the existence of the fast dynamo plasmas in Riemannian space, Phys Plasmas 
15 (2008) in press. 

[14] S. Friedlander, M. Vishik, Chaos 1(2), 198 (1991). 

[15] E. Ott, Chaos in Dynamical systems 2nd Edition (2005) Cambridge. R. Gilmore 
and M. Lefranc, Topology of Chaos (2006), Wiley. 

[16] S.I. Vainshtein, Ya B Zeldovich, Sov Phys Usp 15 (1972), 159. 

[17] G. Paternain, Geodesic Flows Birkhauser (1999). 

[18] C. Chicone and Yu Latushkin, Evolution Semigroups in Dynamical systems and dif- 
ferential equations, American Mathematical Society, AMS-(1999). C. Chicone and 
Yu Latushkin and S. Montgomery-Smith,Comm. Math. Physics 173 379 (1995). 
C. Chicone and Yu Latushkin, Proc of the American Mathematical Society 125, N. 
11,3391 (1997). 

[19] Z. Wang, V. Pariev, C. Barnes and D Barnes, Phys Plasmas 9,5 (2002). 



11 



